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ABSTRACT . 

We have obtained high resolution spectra of V444 Cygni with the International Ultra- 
violet Explorer Satellite. These spectra span both eclipses and include one observation 
at third quadrature. Together with seven archival spectra, they provide reasonably 
complete phase coverage for the system. The variations in the P Cygni profiles of the 
He II and N IV lines imply the existence of a low density region in the WR wind. 
This region occupies a relatively narrow range of orbital phase coinciding with the 
highest terminal velocities observed in C IV. We interpret these data as evidence of an 
interaction region separating the winds of the O-star and Wolf-Rayet star. 
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I. Introduction 


The importance of the Wolf-Rayet binary V444 Cygni = HD 183576 (06-8 III:+ 
WN4 AB; Walborn 1974) is that the system affords an almost unique laboratory for 
studying the structure of a massive radiatively driven stellar wind. In consequence of 
the atmospheric eclipse of the primary by the Wolf-Rayet (WR) secondary, the density 
and velocity structures of the constituent stellar winds are encoded in the variation of 
ultraviolet line profiles with orbital phase. V444 Cygni is the only eclipsing WR + O 
binary which is bright enough to be conveniently studied at high resolution with the 
International Ultraviolet Explorer (IUE). 

Radial velocity studies (Wilson,1940,1842; Keeping, 1947; Munch, 1950) have es- 
tablished the orbital parameters: a = 4OR 0 , Mot = 25A/ 0 , and Mwn = 10 M a . The 
most recent photometric solution (Cherepashchuk and Khaliullin 1973; Cherepashchuk, 
Eaton, and Khaliullin 1984, hereafter CEK), based on analysis of the eclipses at wave- 
lengths between 2460 A and 3.5p, yields the orbital inclination, i = 78°, radius of the 
primary. Rot = 1OR 0 ; and the radius of the WR core, Rw\ = 2.9 R e . These values 
are consistent with earlier photometric determinations (Kron and Gordon 1943, 1950) 
and with a polarimetric determination of the orbital inclination by Rudy and Kemp 
(1978). 

The unequal durations of primary and secondary eclipse have generally been 
considered to imply that the WR core is surrounded by an extended, optically thin, 
strongly limb-darkened outer atmosphere (Kopal, 1944) of radius roughly three times 
that of the core (Russell, 1944). Kopal and Shapley (1946) predicted that electron 
scattering was the principal source of continuous opacity in this dilute envelope, a fact 
confirmed by the observation that the shape of primary eclipse is independent of the 
wavelength of observation (Hiltner, 1949). 


The mass loss rate has been determined by two methods. Bieging, et al. (1982) 
and Abbott, et al. (1986) have measured the free-free emission at 4885 MHz. and 
derived M = 1.4 x 10' 5 M o yr' 1 on the basis of the theory of Wright and Barlow 
(1975). Khaliullin (1974) and Kornilov and Cherepashchuk (1979) have used the rate 
of change of the orbital period (see eq. 1) and the assumption of radial outflow from 
the WR component to infer M = 1.2 x 10“ s A/ 0 yr“ 1 . 

CEK find that the color temperature of the WR core in the optical-ultraviolet 
range is approximately 20,000 K, while the electron temperature is 80,000 - 100,000 
K. The wide wavelength range of their data allowed CEK to determine the run of the 
volume extinction coefficient with radial distance from the WR component. Attributing 
the extinction to Thompson scattering and assuming that the electrons arise entirely 
from the double ionization of helium, they determined that at the WR core, n, = 
9 x 10 ,J cm~ s and v = 400 kms“ 1 . They find that the wind velocity increases steeply, 
attaining an outflow velocity near 2000 kms -1 at r = 1O7? 0 . The wind accelerates more 
gradually at larger radii, attaining its terminal velocity (2500 kins" 1 ) at a distance of 
20 to 30 R 0 . 

The optical spectrum is dominated by high excitation emission lines of He II, 
N III; N IV and N V which display a strikingly diverse range of phase dependent 
variations (Munch, 1950). The blend of He II and N III at A4542 is stronger at primary 
eclipse than at secondary and weakest at elongations, while N III A4640 is stronger 
at secondary minimum than at primary. The flux in the N III A3483 emission feature 
is roughly independent of orbital phase except at primary eclipse, when the emission 
weakens markedly and a strong absorption feature appears. In contrast, N IV A4058 
lacks absorption at any phase and is strongest at primary minimum. 

Kuhi (1968a) obtained narrow-band light curves through secondary eclipse in 
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the principal optical emission lines. His results confirm the diversity of line variations 
(some lines do not show eclipses at all) and also reveal the magnitude of intrinsic in- 
tensity fluctuations in the optical emission lines. Flux differences at constant phase 
of order tens of percent were found, even when the observations were separated by 
only a few orbital periods. Kuhi determined from his observations that the WR atmo- 
sphere, though optically thin in the continuum, is optically thick in the emission lines 
and suggested that electron scattering of line radiation contributes significantly to the 
(fluctuating) characteristics of the narrow band eclipse curves. He was accordingly 
unable to account for his light curves on the basis of geometrical considerations alone. 
From the general tendency for emission lines to be weaker at primary than at secondary 
eclipse, he inferred that a significant fraction of the flux emitted in lines comes from a 
small cylindrical region between the two stars. 

Kuhi (1068b) and Bisiacchi, et al. (1982) found that the profile of He 11 A4686 
is symmetrical at both elongations but asymmetrical to opposite sides during primary 
and secondary eclipse; it is shaded to the blue at primary eclipse and shaded to the red 
at secondary. Sahade (1058) reported the detection of a narrow emission component 
superimposed on the broad A4686 emission feature. He interprets the radial velocity 
variations of this component as indicative of an origin near the Lagrangian point, L t 
of the system. Underhill (1966) has suggested that the presence of sharp displaced 
absorption features particularly near primary eclipse ($ = O.rOO) is evidence for gas 
streaming in the system. 

Recently, low dispersion IUE spectra have been employed to characterize the 
variations of UV emission lines in this system. Eaton, Cherepashchuk and Khaliullin 
(1985, hereafter ECK) noted that the C IV AA1548, 1550 and Si IV AA1396.1403 res- 
onance doublets showed P Cygni profiles with the largest terminal velocities and least 


dependence on orbital phase, while the N IV A1718 subordinate line showed the lowest 
terminal velocity and greatest phase dependence. Additionally, they noted that the He 
II A1640 and N IV A1718 lines showed weak absorption components near secondary 
eclipse. Synthetic photometry in these lines showed: 

1. Primary flux minimum for He II, N IV, and Si IV are displaced from the 
continuum mid-eclipse phase, occurring at t = 0.05. 

2. In He II and N IV, the flux increases between phases 0.25 and 0.45. 

3. The emission in all three of these lines decreases at secondary eclipse. 

They interpret these light curves in terms of a model in which the emitting region 

lies within 9 ± 4 H of the WR component but the absorption arises in a spherical shell 
of inner radius ss 10i? a and outer radius external to the system. The brightening of 
A1640 and A1718 between $ = 0.25 and $ = 0.45 (s= 5 percent of the continuum flux) is 
interpreted as a binary reflection effect: the O star radiation reflected off the expanding 
WR envelope. In contrast, light curves in the N IV) A1486, C IV AA1548.50 and N V 
A1240 display little if any phase dependence, whereby ECK infer that they are formed 
exterior to the formation region of the He II and N IV lines. 

Koenigsberger and Auer (1985) measured fluxes and equivalent widths in emis- 
sion and absorption components of some P Cygni lines as well as velocities of identifiable 
points on the line profiles. They determined that both the emission and. absorption 
components of He II AI640 and N IV A1718 profiles vary with phase. Most striking was 
their observation that the blue-shifted absorption of C IV A1550 is narrower at primary 
eclipse than near secondary eclipse. They attributed this difference to the contribution 
of an O-star wind which reaches a higher terminal velocity than that of the WR wind. 
Accordingly, they suggested the likelihood of a wind-wind interaction region between 
the two stars. 
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Favoring an alternative viewpoint, Underhill (1983a, b, 1984 a,b;1986) has persis- 
tently criticised the widely accepted notion that WR stars possess relatively small hot 
cores and powerful, radiatively driven winds. Underhill and Fahey (1987) argue that, 
because of the presence of weak unresolved emission lines, the UV light curves derived 
(as those of ECK) from low resolution IUE spectra do not represent true continuum 
fluxes. They derive continuum fluxes from a series of 13 high-resolution IUE spectra of 
V444 Cygni. Comparing these data with the optical light curves of CEK, they argue 
for the presence of sufficient phase dependent “third light" shortward of 1500 A to in- 
validate the light curve analysis of Cherepashchuk and his collaborators. On the basis 
of the orbital phase dependent variations of the He II A1640 profile, they recommend 
a model in which the WR component has T,// — 30,000 K, R = 9f?0 and possesses 
large magnetically-supported plumes which contribute much of the optical and UV line 
emission. 

The original goal of the observations reported here was to determine the density 
and velocity structure of the WR wind directly from the absorption components of the 
ultraviolet P Cygni profiles. By subdividing the WR wind absorption profiles through 
secondary eclipse into bins of ~ 300 km s" 1 , we planned to generate a set of light curves, 
each specific to a single projected expansion velocity. We would thus obtain the total 
optical depth along iso-velocity surfaces for each impact parameter. In principle it is 
possible to obtain, from the phase dependence of this optical depth, both the density 
and velocity structure of the WR wind. In practice, the intended Doppler tomography 
was incompletely realized because O-star absorption profile is not separable in the IUE 
spectra, even at quadrature when the velocity difference between the components is 
greatest. Nevertheless, the spectra obtained for this program reveal very important 
features of the binary system. 


II. Observations 

The new observations discussed here were obtained using the SWP camera of the 
IUE satellite during low noise (US1) shifts in 1985 May /June. The large aperture was 
used for maximum throughput and because of the uncertain flux calibration of small 
aperture high resolution spectra. The choice of exposure times was determined by the 
need to achieve maximum phase resolution consistent with acceptable S/N at C IV 
A1550. Because the star is relatively faint and heavily reddened, the longest exposure 
which could be accommodated in a half-shift yielded S/N in the continuum at A1550 
of roughly 5:1, only half that expected for an optimal exposure. These spectra were 
supplemented by those archival images which were taken through the large aperture 
and which have acceptable noise levels. Several archival imageB were rejected because 
of poor S/N, although their omission left phase gaps. A journal of the new observations 
is presented in Table 1, while corresponding data for the archival Bpectra appears in 
Table 2. 

The spectra are phased according to the ephemeris of Kornilov and Cherepashchuk 

(1979) 

H JD(primaryeclipsc) = 2441164.337 + 4. J 212435(£) (1) 

supplemented by a correction for the observed rate of period change (Cherepashchuk 
1982). For SWP 14553, (O-C) = -O.'Oll, while (O - C) = -O.'OIS for all the other 
images. 

The data were analysed using routines available at the Goddard Space Flight 
Center Regional Data Analysis Facility. Because the spectra are underexposed, smooth- 
ing beyond the optimal 3-point standard (which matches detector and spectrograph 
resolution) was required. This is a delicate procedure, as excessive smoothing can not 
only obscure important features but also introduce spurious ones. We performed a 
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number of comparison tests with a running boxcar filter, using the reseaux and the 
profiles of the sharp interstellar C IV lines as fiducial indicators in this connection. 
Easily perceptible changes in profile slope are produced by smoothing to more than 10 
points; serious loss of information begins to occur with 20-point smoothing. We employ 
either 5- or 10-point smoothing, depending upon the S/N of the original image. Those 
archival spectra which were rejected because of inadequate S/N were found to require 
smoothing to 20 or more points before they could be intercompared. 

In the cases of the strongest lines, N V AA1239.43 and C IV AA1548, 51, improper 
interorder subtraction yields negative flux values at the bottoms of the saturated ab- 
sorption troughs. The equivalent widths for these lines were computed by setting all 
negative fluxes to zero. The unsaturated He II A1640 and N IV A1718 absorption lines 
were free of this defect. 

All the lines discussed below exhibit important periodic changes in line structure, 
especially near the cross-over point between the emission and absorption components 
of their P Cygni profiles. In many cases the line profile slope is greatest at precisely 
this cross-over point. These facts make it difficult to define the appropriate continuum 
level which separates the absorption and emission components, especially in relatively 
noisy spectra. The local pseudo-continuumoutside the line is much more easily defined, 
whereby the wavelength interval for the total equivalent width (emission + absorption) 
can be dependably selected. Accordingly we quote only latter quantity; any interpre- 
tation of these data which depends upon the detailed parsing of the P Cygni structure 
is far less certain. 

Finally, we infer the terminal velocity from the wavelength at which the extrap- 
olated profile merges with the local pseudo-continuum. Careful tests of sample profiles 
subjected to different degrees of smoothing confirmed that negligible error in the mea- 


sured terminal velocity is introduced by the smoothing procedure unless the width of 
the smoothing filter significantly exceeds that used in our reductions. 

III. Results 

Our analysis is based upon the variation of the four strongest P Cygni lines 
with orbital phase. The A1640 and A1718 lines, because they arise from highly excited 
states, are more sensitive indicators of wind structure within the binary orbit than are 
the resonance lines. The monotonic increase of SWP camera sensitivity with A > 1500 A 
multiplies this advantage. Conversely, the resonance lines (AA1239.43 and AA1548,51) 
probe the velocity structure at large radii. The measured equivalent widths for the 
principal lines discussed in this section are listed in Table 3. 

a. He II A1640 

The emission line always dominates this profile, whereby the total equivalent 
width is negative. The absorption is strongest at primary eclipse, weak at secondary 
eclipse and completely absent at phase 0.60. The single wave variation of IV^o with 
orbital phase is shown in Figure 1. 

Near primary eclipse the emission profile is rounded, symmetric and red-shifted 
from the rest wavelength in the WR frame. The saturated absorption profile is also 
rounded, mirroring the emission profile shape. Figure 2 shows the contrast between 
the profile at primary eclipse and that at ♦ = 0. p 60. 

Figure 3 shows gallery of line profiles which span the half of the orbit centered 
on secondary eclipse. Near quadratures, the emission profile strengthens slightly, de- 
veloping an overall triangular shape peaked at the WR velocity with a steepened blue 
edge. The absorption profile also becomes triangular, giving a sawtooth appearance to 
the P Cygni line. Near secondary eclipse, the emission profile becomes slightly broader. 
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distinctly plateau-shaped with the steep blue edge characteristic of the quadrature pro- 
files, and symmetric with respect to the WR velocity. These variations of the emission 
line profile do not extend to the high velocity red tail. Longward of approximately 
A1645, the emission profile is invariant with orbital phase. 

The terminal velocity inferred from the absorption line is clearly variable on time 
scales ranging between hours and days. There is no evidence for coherent variation with 
orbital phase, however. The terminal velocities, transformed to the binary center of 
mass frame (-y(06) = 10 km s' 1 ), are shown in Table 4 and plotted against orbital 
phase in Figure 4. 

b. N IV A1718 

A prominent absorption component, visible at most phases on the P Cygni line, 
reaches its maximum at primary eclipse and is weakest near secondary eclipse. The 
emission component is strongest near quadratures, when the absorption is also strong. 
This alternation between absorption- and emission-dominated profiles contrasts sharply 
with the character of the He II line, which is always strongly dominated by emission. 
The total equivalent width of A1718, however, varies in phase with that of He II (see 
Table 3 and Figure 5). 

The N IV emission feature near primary eclipse is round-topped on some orbits, 
flat-topped on others, but always centered significantly redward of the WR velocity. 
When the emission is weakest, the absorption sometimes extends redward of the WR 
velocity. The most important changes in the profile occur near secondary minimum, 
when the absorption component weakens, completely disappearing at 9 = 0. p 52. The 
emission component also weakens beginning at ♦ = 0. p 46, but recovers by 4> = 0. p 76 
The occurrence of this almost pure emission profile is thus not phase-symmetric about 
secondary eclipse. 


Near first quadrature (SWP 24731; ♦ = 0. p 37), the emission and absorption 
profiles are both triangular. The asymmetric emission is peaked at the WR velocity, 
shaded to the red and with a steep blue edge. Near third quadrature, the emission 
profile differs only in that it is shifted redward and truncated, forming a plateau-shaped 
profile. 

The terminal velocity (see Table 4 and Figure 4) is generally consistent with 
that of He II, approximately constant on short timescales at about 2100 km s' 1 . The 
absorption profile is generally asymmetric, displaying the steepest blue edge at the 
quadratures. When strongest the absorption is nearly saturated. An extensive gallery 
of A1718 line profiles was presented in a preliminary report of these observations (Brown 
and Shore 1686). In Figure 6 we show the consistency of the profiles near primary 
minimum and the contrasting secondary eclipse profile. 

c. C IV AA1548.&1 

The total equivalent width of this line iB almost alwayB approximately zero (equal 
fluxes in emission and absorption), suggestive of formation by scattering of light from 
a small core embedded in an optically thick primarily scattering-medium. The deepest 
absorption is seen during primary eclipse. 

The profile variations are not striking at most phases, so that, as noted by 
ECK, the low dispersion spectra supply little information. The emission profile varies 
between two relatively stable shapes. Between phases 0. p 37 and 0. p 60, the emission is 
skewed slightly to the blue. A steep blue emission edge is joined to a correspondingly 
steep red edge of the absoption profile. Between 0. p 62 and 0. p 10, the emission profile 
is shallower and more symmetrical, blending into the absorption line with a relatively 
shallow slope. Figure 7 displays the constancy of the profile through primary eclipse. 


11 


12 



The transition between these forms occurs in a phase interval of approximately 0. p l 
near third quadrature; the phase gap at first quadrature prevents us from observing the 
opposite transition. Except at those phases near secondary eclipse when the absorption 
to emission transition is very steep, the absorption and emission components join at 
1548 ± lA in the WR frame. 

The equivalent width varies by ±4 A (see Figure 8). The mean terminal velocity 
is 2620 ± 420 km s' 1 (3.3o variation). The most dramatic changes occur in the blue 
absorption wing near secondary eclipse. A high velocity blue-shaded tail develops at 
this phase, reaching a terminal velocity of 3800 km s' 1 . The feature is present at phase 
0. p 44 and has disappeared by phase 0. P 7G. Figure 0 shows the phase variation of the 
terminal velocity, while a gallery of line profiles spanning secondary eclipse is shown in 
Figure 10. 

The development of the high velocity tail is accompanied by a change in the 
strength of the emission feature, which becomes slightly broader and shallower. While 
this last effect may be real, it is more likely due to diminished low velocity absorption, 
as the redward edge of the blue-shifted absorption component is steepest at this phase. 
While the profile remains saturated at all phases, the total equivalent width of the line 
becomes most strongly negative (emission dominating) at this phase rather than at 
primary eclipse. 

d. N V AA1239.43 A 

Because this doublet is badly blended with Lyo and falls in a noisy portion of 
the spectrum, the variations can be only roughly characterized. The structure of the 
line is typical for the resonance line of a WR star. The profile is always saturated but 
the structure of the absorption trough, unlike that of C IV, is entirely dominated by 


noise and blends. The terminal velocity appears to be comparable to that exhibited 
by the C IV line, and the same high velocity tail seems to develop at phase 0. p 60. The 
poor S/N of the data in this part of the camera, however, renders these judgments 
uncertain. 

III. Discussion 

The C IV and N V resonance lines certainly contain contributions from both 
components, but the He II and N IV lines appear in the spectra of 06 stars as un- 
saturated absorption lines with no hint of P Cygni emission (Walborn et at., 1986). 
Therefore, all the emission and much of the absorption observed in these lines must 
arise in the WR wind. We therefore discuss the interpretation of the resonance and 
non- resonance lines separately. 

a. He 11 and N IV 

Assuming the CEK solution for the O-star and WR core radii, the well-determined 
orbital inclination implies that only slightly more than half of the WR core is eclipsed 
at secondary minimum. Accordingly, the more than two-fold weakening of the A1718 
and A1640 absorption at that phase implies that either, (a) the CEK value of the WR 
core radius is too large, so that the WR core is actually more substantially eclipsed 
or, (b) these absorption lines are formed almost entirely by absorption and scattering 
of the O-star (rather than the WR) continuum in the WR wind. We adopt the CEK 
photometric solution and hypothesis (b), arriving at the following deductions. 

1. The presence of the absorption features at quadrature implies that the WR 
wind envelopes the entire system. ECK reached a similar conclusion on the basis of 
their continuum light curves for A < 1500 A. 
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2. The disappearance of the A1718 and A 1640 absorption near secondary eclipse 
implies a dramatically lower density of the WR wind in the line of sight at those phases. 
ECK noted the weakening of the absorption features in their low dispersion spectra 
and interpreted it as a result of a binary reflection effect. The reflection (presumably 
resonant backscattering) of O-star continuum light by the (approaching) near side of the 
WR wind was assumed responsible for Ailing in the absorption line and for augmenting 
the emission profile, thus accounting for the observed hump in their narrow band light 
curves. Our high resolution spectra are clearly inconsistent with this suggestion. It is 
well known that the Doppler shift imposed by a "moving mirror" is twice the velocity 
of the mirror. CEK determined that the wind emerges from the WR core at 500 km s' 1 
and reaches terminal velocity (==2000 km s~ 1 in the He II and N IV lines) at roughly 
half the separation of the binary components. Resonant backscattering would thus 
produce a broad enhancement of the continuum, extending between 1000 km s' 1 and 
4000 kins' 1 (4 - 20A) from the rest wavelength in the WR frame. The gallery of A1640 
spectra in Figure 3 show no hint of such a feature. It is difficult to imagine a scattering 
geometry which would so convincingly mimic a reduced column density in the line of 
sight. 

3. The fact that He II and N IV emission is strongest at quadrature implies 
that the strongest emitting/scattering region is on the side of the WR wind facing the 
O-star. The weakness of N IV emission near primary eclipse supports this inference. 

4. The remarkable constancy of the He II emission profile through primary 
eclipse confutes the suggestion of a material stream from the WR toward the O-star, 
as such a stream would presumably produce a recognizable redshifted absorption at 
these phases. The eight profiles which fall in the interval 0. p 89 and 0. p 09 are identical 
within the photometric uncertainty. Redward of 500 km s' 1 from the WR velocity 
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(roughly the velocity which CEK attribute to the WR wind at the photosphere) even 
the profiles at 0. p 76 and 0. p 00 are identical. 

b. CIV 

Notwithstanding the uncertainty in identifying the cross-over pointbetween the 
emission and absorption components of the line profile, it is clear that the red absorption 
edge is never less than 200 km s’ 1 from the WR velocity at any phase. The absence 
of low velocity absorption suggests that the WR wind is optically thick and supersonic 
even at the photosphere, as inferred by CEK from their IUE spectrophotometry. 

The C IV profiles at third quadrature and at d = O p 37 differ very little in 
shape; they are almost identical except for a relative velocity shift of approximately 
400 km s *, somewhat less than the relative WR velocity difference between those two 
phases. It appears that the contribution of the O-star to the P Cygni profile at these 
phases must be slight, whereby the profile variations must be interpreted in terms of a 
non-spherically symmetric wind structure. 

The interval of orbital phase during which the high velocity tail develops in C IV 
and the absorption almost disappears from the A1718 line is moderately well resolved 
in our data. It is evidently asymmetric about secondary eclipse ( 0. p 44 < 4> < 0. p 60) 
and appears to be sharply bounded in phase: one edge occurs between 0. p 37 and 0. p 44, 
the other between 0. p 60 and 0. p 76. The absorption disappears from A1640 not at 
secondary eclipse, but at ♦ = 0. p 60. These characteristics suggest a dynamical origin 
for the phenomenology. 

c. Wind - Wind Interaction 

The simplest explanation is that we are seeing a cavity trailing behind the O 
star: a shadow in the WR wind, produced by the interaction of the winds of the two 

16 



stars. If Walborn’s (1874) luminosity classification (06-8 III) is correct, the O-star 
likely possesses a more substantial wind than typical for a main sequence star of its 
mass. A standing shock would accordingly be expected at the interface between the 
two winds. We interpret the high velocity tail of the C IV resonance line as evidence of 
an oblique shock at the boundary of the interaction region and the absence of the He 
II and N IV absorption as evidence that the O star effectively shadows the WR wind. 
When the cavity crosses our line of sight, we see the lower density O-star wind. In 
the resonance lines, we see the asymptotic material velocity in the interaction region 
projected onto the line of sight. The emission is, as usual, interpreted as the result of 
scattering the O-star light off the WR wind. 

The relative orbital velocity of the O-star with respect to the WR wind yields 
an aberration angle of approximately 0.t)3 between the line of centers and the axis 
of the “wind shadow” cavity. This value is consistent with the observed phase offset 
of the spectral indicators with respect to secondary eclipse. The opening angle of the 
wind cavity appears to be approximately 60° so that the line of sight is parallel to the 
asymptotic conical surface of the interaction region at only two phases: shortly before 
and after secondary eclipse. The resonance lines should exhibit the highest velocity 
blue-shifted absorption at these phases. In the interval O.MS - 0. r 60, the “terminal” 
velocity observed in these lines will be reduced by eosj(90° - i) - 30°] =: 0.9. Figures 
11 and 12 show a pair of cartoons which represent the geometry suggested by these 
considerations. Further details of the model are discussed in the next section. 

Underhill and Fahey (1087, hereafter UF) have proposed a very different model 
for V444 Cygni on the basis of a diverse set of high resolution IUE spectra, many 
obtained during the higher noise US2 shift. An important element of their critique 
of the CEK photometric solution was noted above: the imputed presence of a phase- 
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dependent "third light" in the system. The ultraviolet continuum fluxes which UF 
derive from high resolution spectra are certainly less affected by unresolved weak emis- 
sion lines than are the CEK fluxes. The high resolution calibration, however is far 
less secure than that of the low resolution spectra used by CEK. The high resolution 
photometric calibration is poorest at the short wavelength end of the SWP camera, 
where crowding of the echelle orders complicates the background subtraction. The 
fluxes derived from even optimally exposed spectra are dependable to no better than 
10%. For non-optimal images, other difficulties diminish the photometric accuracy. 
Oliversen (1984) found that near 1300A, where UF find the greatest divergence from 
the CEK light curve, linearity errors of 10 - 20% are not uncommon in trailed low 
resolution spectra if, like the spectra of UF, they are underexposed and degraded by 
high particle background. High resolution spectra may be expected to suffer (at least) 
comparable linearity errors. 

We have supplemented our US1 observations with those of the UF spectra which 
have roughly comparable S/N. Yet, we find no evidence in the resulting set of A1640 
or A1718 profiles for the compact plasma tubes which they suggest link the two stars. 
For example, the peaks of the A1640 profiles in Figure 3 track the WR star velocity. 

One possible origin of the features which they discuss is the smoothing procedure 
which they employ. It is important to note that a 51-point smoothing is not equivalent 
to binning with 10A resolution because the smoothing algorithm performs a running 
average centered on each pixel. In order to improve S/N at the expense of resolution, 
it iB necessary to average the data in non-overlapping wavelength bins. This procedure 
avoids the generation of correlated errors in independent pixels. 

IV. Toward a Model for V444 Cygni 

Interacting winds have been received relatively little attention in the literature. 
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Cherepashchuk ( 1076) and Prilutskii and Usov ( 1076) have presented qualitative dis- 
cussions applicable to WR binaries. Chapman (1081), Huang and Weigert (1082) and 
Girard and Willson (1087) have undertaken more extensive analysis, but for systems 
which bear little resemblence to V444 Cygni. 

For simplicity, we assume that the winds of both stars have reached terminal 
velocity prior to collision. The most recent models of Pauldrach et al. (1086) suggest 
that this is the case for the WR wind in V444 Cygni, but it may be only approximately 
true for the O-star wind. The location of the stagnation point is determined by the 
equality of the wind dynamical pressures. Taking a = (AfjU}/M|V |)' /3 gives x = (jfj) 
for the position of the stagnation point with respect to the O-star in terms of the 
binary separation. Since the velocities of both winds^Pvery high, we neglect the effects 
of orbital motion except for a slight aberration from the line of centers. There is no 
evidence that the WR wind is directly impinging on the surface of the O star, so 
that a must be sufficiently small that the stagnation point is well separated from its 
photosphere. The opening (full cone) angle of the shock relative to the WR star is 
about 6 — 40° so that an approximate solution gives z =s 0.3. 

This places the stagnation point close to the center of mass of the system and 
serves as a constraint on the mass loss rate from the O star. The derived value is 
a = 0.3 — 0.5. Assuming that the wind velocities are approximately equal, Mo = 
(1 - 3) x lO* 6 M 0 yr''. This value is far closer to that of an evolved than main sequence 
O star, consistent with the spectroscopic classification of W'alborn (1974). We also note 
that none of the ultraviolet spectral features conflicts with Walborn’s luminosity class 
assignment. 

IV. Conclusion 

We interpret the abrupt increase near secondary eclipse in the A1550 terminal 


velocity and simultaneous decrease in the material density indicated by A 1640 and 
A1718 absorption as evidence for a shock dominated wind-wind interaction between 
the binary components of V444 Cygni. Our high resolution IUE spectra provide no 
evidence for the existence of material streaming between the two stars, as suggested by 
Kuhi (1968a) or of magnetically configured plasma plumes and tubes, as suggested by 
Underhill and Fahey (1987). 

Because the momentB of the two winds are large and the hydrodynamics of 
their collision is complex, only a preliminary outline of a model is currently possible. 
A qualitative prediction is, however, implied. Of the four other currently known (WR 
+ O) eclipsing binaries, only CQ Cep should exhibit similar phenomenology in IUE 
spectra. The separation of the components of CV Ser is likely too large in comparison 
with the wind momenta for a similar structure to be produced. GP Cep is a relatively 

compact quadruple system in which the winds of all four components may play a role. 

i 

CX Cep is much too faint for IUE observations. Stickland et al. (1984) noted an 
apparent depletion of He II and N IV at primary eclipse, when the leeward side of the 
presumed O-star companion crosses the line of sight. We predict that high resolution 
observations will reveal an abrupt increase in the A1550 terminal velocity at phases 
near the eclipse. 

A second implication of the proposed model is that detectable linear polarization 
would be expected in the UV emission tines because of the presumed importance of 
resonance scattering in their formation. The optical thickness of the lines contrasts with 
the optically thin continuum, rendering a quantitative prediction beyond the scope of 
this paper. 

The importance of V444 Cygni as a unique astronomical laboratory for the study 
of steady state thermodynamics and hydrodynamics in colliding stellar winds amply 
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justifies the theoretical and observational effort which will be required to refine the 
model and test its predictions. 
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Table 1. Journal of New Observations 
SWP | Exposure (minutes) | JD - 2,440,000 [ $ 


25962 

192 

6202.884 

0.098 

25990 

245 

6207.023 

0.080 

25994 

200 

6208.889 

0.523 

26000 

192 

6209.891 

0.761 

26007 

200 

6210.886 

0.997 

1 26008 

205 

6211.248 

0.036 

26030 

200 

6214.890 

0.948 

26031 

205 

6215.051 

0.986 

26041 

202 

6217.050 

0.461 

26065 

200 

6218.812 

0.879 

26066 

200 

6218.968 

0.916 
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Table 2. Journal of Archival Spectra 


SWP 

Exposure (minutes) 

JD - 2,440,000 

4> 

14553“ 

90 

4811.396 

0.773 

24731* 

170 

6056.580 

0.366 

24734* 

180 

6057.578 

0.603 

24740* 

195 

6058.587 

0.843 

24745* 

200 

6059.583 

0.079 

25904' 

200 

6195.887 

0.437 

25905' 

205 

6196.049 

0.475 


Notes: 

• VJLSPA image. 

* US2 image. 

' US1 image. 
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Table 3. Equivalent Widths (A) Table 4. Terminal Velocities (km s' 1 ) 


SWP 

« 

N V 1240 

C IV 1550 

He 11 1640 

N IV 1718 


SWP 

* 

IV ( 06 )* 

tv ( W , R )* 

v ,( 1550 ) 

Vi ( 1640 ) 

Vi ( 1718 ) 

14553 

0.773 

2.8 

- 2.1 

- 23.9 

- 2.4 


14553 

0.773 

110 

-292 

-2380 

-2040 

-2300 

24731 

0.366 

0.48 

- 0.8 

- 18.9 

- 3.7 


24731 

0.366 

-90 

238 

-2340 

-2050 

-1590 

24734 

0.603 

- 3.0 

2.8 

- 25.4 

- 5.8 


24734 

0.603 

72 

-174 

-3480 

-1170 

-1810 

24740 

0.843 

2.7 

- 1.8 

- 20.2 

0.8 


24740 

0.843 

100 

-244 

-2270 

-1890 

-2570 

24745 

0.079 

3.0 

0.6 

- 15.7 

3.1 


24745 

0.079 

-58 

155 

-2100 

-1750 

-2520 

25904 

0.437 

1.2 

- 0.5 

- 19.2 

- 3.3 


25904 

0.437 

-46 

128 

-3230 

-1780 

-1930 

25905 

0.475 

4.2 

1.1 

- 17.1 

- 3.5 


25905 

0.475 

-19 

58 

-3050 

-2330 

-1640 

25962 

0.098 

1.3 

1.8 

- 16.1 

2.1 


25962 

0.098 

-69 

186 

-2180 

-2400 

-2100 

25990 

0.080 

2.8 

3.1 

- 15.3 

3.3 


25990 

0.080 

-58 

157 

•2180 

-2170 

-1820 

25994 

0.523 

5.3 

3.5 

- 20.7 

- 4.4 


25994 

0.523 

17 

-34 

-3370 

-1980 

-2180 

26000 

0.761 

4.3 

- 3.3 

- 21.3 

- 3.2 


26000 

0.781 

120 

^,< J67 ^ 


„< tfS 60 ) 


26007 

0.997 

- 3.7 

0.7 

- 19.3 

1.7 


26007 

0.997 

2 

4 

-2410 

- 1910 " 

-1910 

26008 

0.036 

- 12 . 2 : 

- 1.0 

- 14.3 

3,0 


26008 

0.036 

-27 

78 

-2310 

-1880 

-1980 

26030 

0.948 

- 21 . 1 : 

2.6 

- 12.8 

4.7 


26030 

0.948 

38 

-88 

-3290 

-2760 

-2560 

26031 

0.986 

- 6.8 

0.6 

- 15.0 

1.8 


26031 

0.986 

10 

-17 

-2290 

-2440 

-2530 

26041 

0.461 

1.7 

- 1.1 

- 21.1 

- 2.2 


26041 

0.461 

-29 

84 

-3810 

-1980 

-2030 

26065 

0.879 

- 3.7 

- 0.0 

- 13.9 

1.0 


26065 

0.879 

83 

-200 

-2310 

-2440 

-2480 

26066 

0.916 

- 6.5 

0.6 

- 13.1 

4.0 

. 

26066 

0.916 

60 

-144 

-2260 

-2340 

-2250 


Note: 

“ Derived from Munch (1950). 
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FIGURE CAPTIONS 


DOUGLAS N. BROWN:' Department of Astronomy FM-20, University of Washington, 
Seattle, WA 98195 

STEVEN N. SHORE: Astrophysics Research Center, Department of Physics, New Mex- 
ico Institute of Mining and Technology, Socorro, NM 87801 


figure 1 . 


Figure 2. 


Figure 3. 


Figure <i. 

Figure 5. 
Figure 6. 


Figure 7. 

Figure 8. 
Figure 9. 

Figure 10. 


Figure 11 . 


Figure 12. 


Orbital phase variation of He II Xi690 total equivalent width. 

The uncertainty in phase is less than the size of the plotting 
symbols. 

Comparison of four He II X.1690 line profiles near primary 
eclipse with that at 0^60 (SWP 29739). The eclipse profiles 
are from SWP 26007, 26008, 26030 and 26031. Note the complete 
absence of absorption at 0?60. 

Gallery of He II X.169Q line profiles spanning secondary eclipse. 
Phase advances from the bottom (0?37) to the top (0?89) of the 
figure. Note that the absorption component Is almost absent between 
o!99 and 0?76. 

Orbital phase variation of terminal velocities measured In 
He II X1690 (trloLv.jI'i and N IV Xl718 (dots). The uncertainty 
in phase is less than the size of the plotting symbols. 

Orbital phase variation of N IV X1718 total equivalent width. 

(a) Overplot of N IV "Xl 718 line profiles near primary eclipse. 

The spectra are taken from the same images as In Figure 2. 

(b) Comparison of N IV Xl718 line profile at secondary eclipse 
ISWP 2599 1 * ■ 0?52) and lust before onset of primary eclipse 
(SWP 26065, 0?88) . Compare with Figure 3; note that the 
absorption component of Xl 71 S disappears at earlier phase 
than does that of X1690. 

Overplot of C IV X1550 line profiles near primary eclipse. 

The spectra are taken from the same Images as In Figure 2. 

Orbital phase variation of C IV Xl550 total equivalent width. 

Orbital phase variation of C IV X1550 terminal velocity. Note 
that the increased wind velocity coincides in phase with the 
weakening of He II and N IV absorption. 

Gallery of C IV\l 550 line profiles spanning secondary 

eclipse. Phase advances from the bottom (0f37) to the top 

(0?84 ) of the figure. Note the presence of enhanced 

blueshlfted absoptlon between oC99 (SWP 25909) and of60 (SWP 26000). 

Cartoon of the proposed wind - wind Interaction in V999 Cygnl. 

Top view along the normal to the orbital plane. The SWP image 
numbers of spectra discussed In the text are Indicated at the 
corresponding orbital phases on the reference circle. 

Cartoon side views of the V999 Cygnl system at secondary eclipse 
showing the proposed structure of the wind • wind Interaction. 
Direction of observer Is to the right. 

(a) View along the line of nodes. 

(b) View along the line of sight. 
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